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Abstract 

We explore CP violating aspects in the Higgs sector of models where new vectorlike quarks carry 
Yukawa couplings mainly to the third generation quarks of the Standard Model. We point out that 
in the simplest model, Higgs CP violating interactions only exist in the hWW channel. At low 
energy, we find that rare B decays can place similarly strong constraints as those from electric dipole 
moments on the source of CP violation. These observations offer a new handle to discriminate from 
other Higgs CP violating scenarios such as scalar sector extensions of the Standard Model, and 
imply an interesting future interplay among limits from different experiments. 
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I. INTRODUCTION 


Now that a Higgs boson has been observed with properties similar to those predicted by 
the Standard Model (SM), the next critical task is a program of precision measurements 
of its properties. Studies of the Higgs mass, coupling strengths, and production and decay 
channels are well advanced. However, less attention has been paid to the possibility of 
observing CP violation in the Higgs sector and the purpose of the present work is to explore 
such a possibility. 

If the 125 GeV boson is measured to be a mixture of CP even and odd states, it imme¬ 
diately indicates that there must be new physics not far above the electroweak scale. One 
of the most straightforward ways is to extend the scalar sector of the SM, and the simplest 
case is the complex version of the 2 Higgs doublet model (C2HDM) [1-6]. The presence 
of CP violation leads to observable changes in Higgs production and decay rates, as well 
as contributions to low energy observables such as electric dipole moments [7]. It has also 
been pointed out that some of the heavy scalar decay channels could be sensitive to CP 
violation [8], and could be probed at the LHC and future colliders. The complementarity of 
LHC and low energy measurements for constraining Higgs CP violation has been explored 
in Refs. [7-10]. 

Alternatively, it is also possible to extend the fermion sector of the SM while keeping 
the scalar sector minimal. The simplest case is to introduce a chiral fourth generation or 
mirror family, but they are strongly disfavored after the discovery of the Higgs boson, since 
they would lead to a large enhancement in the Higgs production rate. The next simplest 
extension is to introduce vectorlikc quarks (VLQs). Vectorlike fermions are defined as having 
the same gauge quantum numbers for both left- and right-handed fermion pairs, and thus 
they do not generate chiral anomalies and their effects decouple in Higgs physics. They 
are the ingredients of beyond the SM frameworks like the little Higgs [11] or composite 
Higgs models [12, 13], and theories of extra dimensions [14-16] or extended supersymmetry 
[17, 18]. VLQs have also been discussed recently in light of modified Higgs couplings such 
as that to two photons [19, 20]. 1 The current LHC lower limit on the masses of VLQs from 
direct searches is around 800 GeV [24] regardless of the decay modes. The indirect effects of 
VLQs in electroweak precision measurements and flavor physics have also been extensively 

1 There are also recent studies which extend both scalar and fermion sectors of the SM [21-23]. 
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explored in the literature [25-29]. 

In this work, we consider the CP violating aspects of the VLQ models, which have been 
less studied. Our goal is to examine their impact on the CP nature and interactions of the 
Higgs boson, as well as the low energy constraints [30]. For simplicity, we focus on the simple 
cases where VLQs are in a single representation of SU(3) C x SU(2)l x U(1)y- We consider 
only fermion representations that can have new Yukawa couplings with the SM quarks and 
Higgs doublet. These new Yukawa couplings could provide a new source of CP violation at 
the weak scale. Under these assumptions, we find that only the case where the VLQ lies in 
the (3,2,1/3) representation can generate significant CP violation in Higgs physics. 

We study CP violating phenomenology in the doublet VLQ model using an effective the¬ 
ory language where the heavy VLQs have been integrated out. Then the CP violating effects 
manifest themselves through a new right-handed charged-current interaction mediated by 
the fU-boson. We clarify the source of CP violation in this model in Section II. In Sec¬ 
tion III, we calculate the loop induced CP violating Higgs interactions with gauge bosons. 
Interestingly, we find that CP violation only exists in the hWW coupling, but not in the 
hZZ[ 31], hj'y, or hZ 7 [32] ones. In Section IV, we explore the current constraints on this 
coupling from low energy measurements, including EDMs and B physics, and comment on 
the future prospects. 

II. CP VIOLATION FROM VECTORLIKE QUARKS 

The vectorlike quark representations that allow new Yukawa couplings with SM quarks 
are summarized in Table I 2 . 

The key point we want to make is that only the doublet (T, B) model can offer non-zero 
(unsuppressed) CP violation for the Higgs boson, through the Yukawa coupling to third 
generation quarks 3 . To see this, we write the Yukawa sector of this model, 

= UtQzLHt R + UbQ3LHb R + MQ' L Q' R + M'Q 3L Q' R + \ t Q’ L Ht R + X b Q' L Hb R + h.c. ,(1) 
where H = (</>+, 0 O ) T is the SM Higgs doublet, H = icroH*, Q 3L = (t^&Q is the third gen¬ 
eration left-handed quark doublet and Q' RR = (T,B) j j R are the vectorlikc quark doublets. 

2 We find that introducing any number of singlet T LR and/or B LR fields does not give rise to new CP 

violating phases since any potential new phases can always be rotated away by a held redefinition. 

3 We assume only one representation of VLQs is present. 
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VLQ models 

Representation 

CP violation 

Tl,r 

(3,1,4/3) 

no 

Bl,r 

(3,1,-2/3) 

no 

C T,B) LtR 

(3, 2,1/3) 

yes 

(X,T) l ,r 

(3, 2, 7/3) 

no 

C b,y) l , r 

(3, 2,-5/3) 

no 

(X, T, B) L)R 

(3, 3,4/3) 

no 

(T,B,Y) l ,r 

(3, 3,-2/3) 

no 


TABLE I. Models of vectorlike quarks and their representations under SU( 3) c x SU(2)r x U(l)y, 
together with the possibilities of introducing new physical CP violating phases. 


In general, the SM gauge invariance permits us to generalize the fields above, (Qsl, A, bn), 
to linear combinations of all three generations. In our study, we assume the fields of Eq. (1) 
are dominantly composed of third generation fermions, because they have the largest Yukawa 
couplings and thus have the strongest impact on CP violation in Higgs physics, which is the 
motivation of this work. To be concrete, we can take advantage of the hierarchical structure 
of the CKM matrix, and define Eq. (1) in the basis where the SM 3x3 blocks of the up- and 
down-type Yukawa matrices are close to diagonal up to CKM-like rotations 4 . This helps to 
suppress the mixing between heavy VLQs and the first two generation quarks and minimize 
the low energy flavor changing effects in the spirit of next-to-minimal flavor violation [33]. 

From now on, we will focus on the mixing between VLQs and the third generation quarks. 
Since Q^l and Q' L have the same quantum numbers, one can always redefine holds and set 
the parameter M' = 0. After electroweak symmetry breaking, the quark mass matrices take 
the form, 


A™ — (A, T l ) 



+ ( b L , B l ) 


5§ o 
vf M 



( 2 ) 


In general the parameters are all complex, and one can remove unphysical phases by re¬ 
defining the phases of the fields. Under the gauge invariant transformations, Q' R —y Q' R e ia , 


Alternatively, we could work in the basis where the 3x3 block of down quark Yukawa couplings is already 


diagonal. With this assumption, the down quark sector is free from new contributions to flavor violating 


processes. 
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Q'l Q'l^i Qsl —*■ Qzlz 1 " 1 i tR —>■ t R e tS , b R —> bRe ia , the parameters change to 

Vt yte i[6 - l] , A t -+ \ t S s ~® , 
y b -+ y b e i{a ~^ , A fe \ b e^~V , 

M ->■ Me i{a ~ p) . (3) 

Clearly, what is invariant is the combination of phases of the parameters arg (y b ) + arg(A t ) — 
arg(yt) — arg(Afc), or the quantity, 

Im(y b \ t yt K) = \ytyb^b\\e ie . (4) 


This is the only new source of CP violation in this model that can enter into physical 
processes. It is convenient to use the above rephasing freedom to rotate the phase 6 into X t 
and make the other parameters real. In this case, from the quark mass terms, y q (v+h)qq/\/ 2, 
we can first assign the Higgs boson to be CP even. Then any coupling between h and CP 
odd operators induced by Eq. (4) indicates CP is violated in the Higgs sector. 

We diagonalize the mass matrices in Eq. (2), and obtain the mass eigenstates 


t R = cos OrPr + sin d^e l6 T R . 


b R = cos 0 b R b R + sin 9 b R B 


R n Ri 


(5) 


where T Rl B R are orthogonal to t R , bR, respectively. There are similar mixings among the 
left-handed fields, parametrized by angles 9 f L and 9 b L . The mixing angles among the right- 
handed quarks satisfy, 


tan 29’ n = — „ ., „ ON . . 

M*-(\? + y?)vy2 


M 


(i = t,b ), 


( 6 ) 


and the last step keeps only the leading term in the v/M expansion. The angle 9 b R denotes 
the mixing between the SU(2) L singlet b R and doublet B R , and is constrained by electroweak 
precision measurements such as Z —y bb [25, 34, 35]. The mixings among left-handed quark 
fields only appear at order (v/M) 2 and are much smaller [25, 26]. 

The lower bound on the VLQ mass scale is around 800 GeV from direct searches at the 
LHC [24], which suggests that we can integrate them out when studying Higgs physics. 
Since T R , B R lie in an SU(2) L doublet, integrating out the heavy vectorlike quarks yields an 
anomalous Wtb interaction, 


£eff — or [ —t= 




(7) 
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where 


■ /i t • /}/) if) ft f I ^ jf) 

a R = sm sm 9 b R e lt> ~ 2M ' 2 e* . (8) 

As discussed above, in this model CP violation must appear in physical processes through 
the combination of couplings, lm(y b X t y^ XI) . The new right-handed Wtb coupling a R ob¬ 
tained here is proportional to AtA£. Therefore, a physical process that makes the CP vio¬ 
lation manifest should involve both left- and right-handed currents in order to allow mass 
(Yukawa coupling yt,y b ) insertions. We write the most general renormalizable Wtb coupling 
as 


£ e // — lPl + a R P n )bW+ + h.c.. 


( 9 ) 


We neglect the hat symbol for mass eigenstates hereafter. In the SM, = Vt b — 1 and 
a R = 0. For the vectorlike quark doublet model we consider, a R is given by Eq. (8), and the 
deviation of a R from 1 occurs only at higher order in v/M. 

Following a similar reasoning, we have also examined other representations of VLQs. 
Interestingly, none of them can offer an irreducible CP violating phase such as that in 
Eq. (4), i.e., under the same assumptions, the Higgs boson is essentially CP even in those 
models. This observation places the VLQ doublet (T, B ) model in a unique place in the 
perspective of CP violation. In the coming sections, we will explore the ( T,B) model as a 
theory of Higgs CP violation, and study in detail its predictions in phenomenology and the 
constraints on the parameters of the model. 


III. CP VIOLATION IN HIGGS BOSON INTERACTIONS 

The Yukawa interactions between vectorlike and SM quarks introduce a new source of 
CP violation. One of the consequences is that the Higgs boson will obtain CP violating 
interactions with the other SM particles. As discussed in Eq. (4), in the model with a single 
VLQ doublet, ( T, B ), the physical CP violating phase has to appear via the combination 
of parameters, \m(y b XtylX b ). This means that the diagram giving CP violating interactions 
to the Higgs boson must involve both top and bottom quarks. As a result, the leading CP 
violation in this model resides only in the hW+ v W~ tlv operator, generated at loop level as 
shown in the left diagram of Fig. 1. At this order, CP violating tree level Higgs-quark or 
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FIG. 1. Feynman diagrams generating CP violating Higgs couplings. The label “L (R)” means a 
left- (right-) handed current interaction with the VF-boson. The left diagram is in the full theory, 
and the right one is in the effective theory when the vectorlike quarks are integrated out. The 
right-handed current Wtb vertex is derived in Eq. (9). For the hWW coupling, one of the 4>° fields 
is replaced by the electroweak vev, and the other replaced by h. 


loop level Higgs-Z-boson and Higgs-photon interactions are absent. The direct probes of 
CP violating hWW interactions at the LHC have been discussed in the h —>■ WW decay, 
the hW associated production channel and the WW fusion channel [36, 37]. 

In the heavy fermion limit, the gauge invariant operator generating the Higgs-gauge boson 
CP violation starts from dimension 8 in this model, 

A = § [e tj H i ( a ‘y„H k WZ „] 


i H m (a b ) n l H l W bfU ' 


( 10 ) 


a a are the Pauli matrices, W^ u = \e^ va pW al3 and i,j,k,l,m,n = 1,2 are SU(2)l indices. 
After electroweak symmetry breaking, H 1 = (f> + = 0 and H 2 = (f>° = (v + h)/y/2 in the 
unitary gauge. This projects out the CP violating hW+ v W~ tlv interaction 


£« ->• hW“,W, a * = cff-WtW; 


A 4 


/W flU 


LLV flU 


( 11 ) 


In the last step, we define the coefficient a^ in the same notation as Eq. (1.15) in the Higgs 
Working Group Snowmass report [38]. 

In reality, because the top quark mass is comparable to the center-of-mass energy a/s in 
the Higgs production and decay processes, the coefficient aY becomes a form factor. We first 
calculate the form factor for CP violating W ( qi)h(p ) associated production [39, 40] via an off- 
shell W*(q 2 ). In this case, the momenta satisfy qf = A7,|- , p 2 = m 2 h and s = q 2 > (m^+M^y) 2 . 
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FIG. 2. Real part of the form factor a^ as a function of \J in units of 10~ 3 x Im(a£a#). For 
Wh associated production, the kinematics require > rrih + Mw, he., the white region to 
the right of the shaded region. For h —» WW* associated production, the kinematics require 
0 < < rrih — M\y , i.e., the white region to the left of the shaded region. 


The leading contribution to the form factor aY can be conveniently calculated in the effective 
theory when the vectorlike quarks are integrated out, as shown on the righthand side of 
Fig. 1. We find 


,W( 


2 \/2N c G F m t m b * 

'<&) ~ - 4 ^ 2 - MWr) 


f*l — X 


dx 


dy{ 1 - x-y ) 


1 


1 

A b 


( 12 ) 


where the denominators are 


A t = (x + y)z 2 t - xyzl + (x + y - l)(xzf + yz\ - zl) , 

A b = (x + y)zl - xyz 2 h + (x + y - l)(a;z? + yz\ - z, 2 t ) , (13) 

and q 2 are the off-shell momenta of the W gauge bosons, z a = m a /M\y, (a = t,h,b), 

z\ = qf/M'w = 1, and z\ = q 2 /M ^ = s/M^. The 1/A^b terms correspond to the diagrams 

where the Higgs field is attached to the top (bottom) quark propagators. There is a minus 
sign between the two pieces in the integrand of the Feynman parameter integral. From the 
analysis of [37] , only the real part of the form factor a^ contributes to the final CP violation 
observable, he., a phase shift in azimuthal angle. In Fig. 2, we plot the real part of a^ as a 
function of ^Jq f. For Wh associated production, the kinematics require > m/, + Mw, 
i.e., the physical region is the white region to the right of the shaded region in the plot. 
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We next examine the form factor in the decay process h(p) —>■ W{q\)W*{cfe). In this 
case, the momenta satisfy p 2 = m\, q\ = M^, 0 < q% < ( rrih — M w ) 2 . The integral of 
the 1/A t term is real. On the other hand, the integral of the 1/A& term has an imaginary 
(absorptive) part, which is due to a pole in y corresponding to the on-shell cut of the bb 
propagators. Numerically, we find the integral over 1/A 4 and dispersive part of the 1/A;, 
integral almost cancel each other. The integral is dominated by the absorptive part of the 
I/A 5 integral, which we find to be of order 1 for all values of q\. Physically, it indicates that 
the CP violating effect in the h —>■ WW* decay is dominated by processes where the Higgs 
boson first decays to bb and then the bb re-scatter into WW* 5 . 

The coefficient a ^ calculated above is proportional to the quantity Im(a](a/j) and is of 
the order ( v/M ) 2 . There is another contribution obtained by changing the right-handed 
current to a left-handed one in the heavy quark vertex (left diagram of Fig. 1), and we have 
checked that this contribution is 0(v/M ) 4 and is subdominant. 

Using the central values of masses and constants from the PDG [41], we find in both 
processes the coefficient for the CP violating hWW interaction is 

aY — 1CT 3 x Im (a* L apt) ■ (14) 

The first factor contains the usual loop factor and the bottom quark Yukawa coupling. The 
CP violating parameter lm(a* L aR) depends on the model parameters AA t . The goal of the 
next section is to explore the current and future experimental constraints (sensitivities) to 
Im (a^cqj). 


IV. CONSTRAINTS 

In this section, we explore phenomenological constraints on the parameter Im(a^aR) 
relevant for the CP violating hWW coupling. We fold the most relevant limits come from 
the electric dipole moments and the rate and CP asymmetry of the rare B decay b —> s'y. 

5 In general a cut is not necessary for CP violation to occur because the final state W + W~ is already an 
eigenstate under CP. The cancelation between the 1/A t and dispersive part of the 1/Ab integrals seems 
accidental. 
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A. Electric Dipole Moments 


Electric dipole moments are sensitive probes of new sources of CP violation. We first 
study the constraint from the electron EDM. The interactions Eq. (9) can contribute to 
the electron EDM through the two-loop Barr-Zee type diagrams as shown in Fig. 3. This 
contribution has been calculated analytically in Ref. [42], 


d e 


a 


8n 2 sin 4 9 W M W 

f 1—xi 


z e z t z b lm(a* L a R ) 


Qb 

x — / ax i 


dxo 


x(x — 1 ) 


log 


9b 


Xg b -x(l-x)f x(l-x) g b -x( l-x)_ 


(b ^ t ), 
(15) 


where x = x 1 + x 2 , z a = m a /M w , (a = e, t, b), and g t = x(z 2 - z b ) + z 2 , g b = x(z 2 - z 2 ) + z\. 
Here d e is the coefficient of the effective EDM operator, 

C e ff D d f , d(7 ll ^f 5 eF ,lu . (16) 

Numerically, we find, 

d e ~ —1.58 x 10~ 2 'Im(a2ai?) e • cm. (17) 


The current experimental upper limit on the electron EDM is |d e | < 8.7 x 10 29 e ■ cm at 
90% CL from the ACME experiment in 2013 [43]. This translates into the upper bound, 


|Im(a/j)| < 0.055. 


(18) 


It turns out that the electron EDM constraint is weaker than the one from B physics, as 
will be discussed in the next subsection, although the EDM constraint will become relevant 
if the current ACME limit is improved by only a factor of a few. In Fig. 5, the horizontal 
magenta dotted line shows the future exclusion if the limit reaches 10 times the ACME-2013 

limit. 

Next, we consider the EDMs of the neutron, the proton and the mercury atom. These 
constraints usually involve large hadronic and nuclear physics uncertainties. However, given 
the future prospects of these experiments, they could become relevant. There are several 
contributions to these observables. The first includes light quark EDMs, from a similar 
diagram as Fig. 3, with the replacement (e, u) —> ( u , d) or (d, u ). At ji = 1 GeV, 

d u (g) ~ —2.3 x 10~ 26 ?/ 1 Im(a^a^) e • cm, 
dd(g) — —4.6 x 10 - 26 ?/iIm(a](a^) e • cm. 
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(19) 

( 20 ) 










FIG. 3. Barr-Zee diagrams that contribute to electron EDM. The crossed diagrams with L e-)- R 
are not shown. 


Here the renormalization group (RG) running effect from the Mw scale down to the GeV 
scale is taken into account, 


Vi = 


a s (M w ) 

16 

23 

a s (m b ) 

16 

25 

a s {m c ) 

_ a s (m b ) _ 


_a s (m c )_ 


_a s (l GeV)_ 


16 

27 


~ 0.417. 


( 21 ) 


Hereafter we have used the NLO values of a s at various scales in the following table. 


a s (Mw) 

a s (m b ) 

a s {m c ) 

a s (l GeV) 

0.120808 

0.218894 

0.382156 

0.455862 


The contribution of Eq. (19), (20) to the neutron EDM, d n ~ —0.35 d u (fi) + 1.4d d (/i), 
is too small to yield a competitive constraint on lm(a* L an) in view of the current limit 
\d n \ < 2.9 x 10 -26 e ■ cm. 

There is no light quark chromo-EDM at one loop level in the VLQ model. Instead, there 
is a contribution to the chromo-EDM of the bottom quark, shown in Fig. 4 6 . The bottom 
quark chromo-EDM can contribute to the EDMs via matching to the three-gluon Weinberg 
operator at a low scale. The effective Lagrangian for the two operators takes the form [44], 

pabc ^p,i/p(T /^a ib \s~ic 
\ Gt,, Ur 




C& 


Ceff D if 2 g s m b ba^T a bG a ^ + *G c pa . (22) 

The calculation of the coefficient of the chromo-EDM operator is similar to that in the 


left-right symmetric model [45]. At the weak scale, 

Sb(Mw) y/2Gp iTit 


A 2 


8n 2 m b 


Im (a* L a R )f(z t ), 


(23) 


6 Unlike Ref. [42], we find that diagrams similar to Fig. 3 but with photon lines replaced by gluon ones and 
leptons replaced by light quarks vanish and do not give rise to Chromo-EDMs. 
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FIG. 4. One-loop contribution to the bottom quark chronro-EDM. 

where f(z t ) = [l — | z\ — \z f + \zl log zj] /(I — z%) 3 ~ 0.35. Interestingly, there is an 
enhancement factor ( mt/m b ) [46]. At the scale m b , the matching condition is C ( ~,(rri b ) = 
j^a s (m b )8 b (m b ) [7]. Taking into account the RG running, the coefficient Cq at /j = 1 GeV 
is, 


CM 

Oi s {m b ) 

a s (M w ) 

14 

23 

Oi s (m b ) 

29 

25 

a s {m c ) 

5 b (M w ) 

4.5 x 10- 9 Im(a£a R ) 

A 2 

12tt 

_ a s (m b ) _ 


_a s {m c )_ 


_a s (l GeV)_ 

A 2 _ 

GeV 2 


(24) 


The final contribution to the neutron EDM is dominated by the Weinberg operator [7, 44] , 

d n ~ (2 x 1(T 20 e • cm) ^ Cq{[i) ~ —5.5 x 10 _ 24 Im(a2aj?) e • cm , (25) 

where we have used the hadronic matrix element given in Ref. [7] . The current limit on the 
neutron EDM, \d n \ < 2.9 x 10 ~ 26 e ■ cm at 95% CL, translates into the upper bound, 

|Im(a*a jR )| < 0.5 x 10 ~ 2 . (26) 

Because the three-gluon operator is an isospin singlet, the proton EDM in this model is equal 
to the neutron EDM. A possible future experiment measuring the proton EDM is expected 
to give a strong constraint [47, 48]. 

The EDM of mercury 199 Hg is also sensitive to the three-gluon operator, which contributes 
through the Schiff moment [7, 44], 

dn g = K S 29 f lN (ao7(o) + «i7(i)) Cq - 3 - 9 x 10 _27 Im (a* L a R ) e ■ cm , (27) 

where we have used the conventions and values of parameters given in Ref. [7]. The current 
limit on the mercury EDM, cffig < 3.1 x 1CT 29 e • cm, translates into the upper bound, 

|Im(a]ia R )| < 0.8 x 10 ' 2 . (28) 
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B. B physics 


There are strong constraints on the parameter I m(a* L a R ) from the b —> s 7 channel, both 
from the total rate and the CP asymmetry. The effective Lagrangian relevant for this process 
takes the form, 


4W = + c a 3 £±s L ^n R G;„ . (29) 

The contribution of new right-handed current interaction to the Wilson coefficients at the 
scale Mw are [49], 


Ac 7 {M w ) — a R -Af 7 (m 2 /M^ v ), 
Ac 8 (M w ) = a R —f 8 (ml/M , 

Tflb 


(30) 


where the form factors are, 


, , , -3a; 2 + 2x , -5a; 2 + 31a; - 20 

M x ) = “777—ttt logx + 


2(x — l) 3 


12(a; — l) 2 


.. , , 3a; , x 2 + x + 4 

hix) =2(^W gX ~ 4(x — l) 2 ' 

When we take into account the 1-loop RG running corrections from M w to /a 
effective coefficients are [50], 


(31) 
mb, the 


c 7 (/i) = a.cfV) + vl 6/23 Ac 7 (M w ) + ^(r/ 6 U/23 - rjl 6/23 )Ac 8 (M w ), 

Cs(fj) = a LC§ M (/i) + J]l 4/23 Ac s (M w ), (32) 


with Cy M (/i) = —0.31, cf M (/i) = —0.15 and rjb = a s (Mw) / a s (mb) 
The B —» X s7 decay rate in the VLQ model is then given by 

2 


B(B -)• X sl ) = B{B -7 W 7 ) sm 


c 7 


r SM 

c 7 




0.55. 


(33) 


The SM prediction has a central value B(B —* A s7 )sm = 3.15 x 10” 4 . The world average of 
the measurements is B(B —» AQ 7 ) = (3.55 ± 0.24 ± 0.09) x 10~ 4 [51]. 

In order to present this result as a limit, we note the fact that in the VLQ model the 
contribution to a R dominates over the deviation a R — 1. Using results from the EDM 
discussions that lm(a* L aR ) is already constrained to be less than ~ 10~ 2 , we conclude that 


lm(a* L a R ) ~ Im(a fl ). 


(34) 
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FIG. 5. Constraints on the Ke(ap) and Im(afl) parameter space from EDMs and the B —> X s 7 
decay. The red shaded region is consistent with the branching ratio B(B —> Xg'y). The regions 
inside the solid yellow circles are excluded by the direct CP asymmetry, Acp, in B —» Xg'y decay. 
The EDM constraints are shown by the blue dot-dashed lines (current neutron EDM, central value), 
green dashed lines (current mercury EDM) and magenta dotted lines (future electron EDM with 
10 times better limit than the current one from ACME). The exclusion regions are in the direction 
of arrows. The neutron and mercury EDM typically involve large nuclear/hadronic uncertainties, 
and here the lines only show the central value of their bounds. 


In this case, the constraint from the B —> X s ~f decay rate measurement is presented in the 
parameter space of Re(a^) versus Im(a/j), as shown by the red shaded region in Fig. 5. 

The direct CP asymmetry in B —s- Xg'y decay rate is [52] 


A cp 


a s (m b ) 




40 A c 
- 1 

9 m b 


(1 + e s ) 



(35) 


where c 2 = 1.11, e s = —0.007 + 0.018i, A c = 0.38 GeV, and a s (rn b ) = a s (M w )/r] = 0 . 21 . 
The most stringent experimental measurement is from BaBar, Aqp = (1.7±1.9±1.0)% [53]. 
Again, we show this as a constraint in the Re(a^) versus Im(a^) plane in Fig. 5. The regions 
inside the yellow circles are excluded. 

Summarizing the EDM and the B physics constraints, we conclude that in the doublet 
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VLQ model it is still possible to have Irn(a^) as large as order 0.01. From Eq. (14), this 
implies the CP violating hWW coupling a}'} is currently constrained to be at most 10~ 5 . 
The next generation EDM search is expected to further narrow down the allowed window 
of Irn(aft). In the case of discovery, this would trigger an exciting interplay between the 
studies of CP violation in a future B factory and a future Higgs factory. 


V. CONCLUSION 

In this work, we have studied the possibility of introducing CP violating interactions to 
the 125 GeV Higgs boson by extending the fermion sector of the SM with vectorlike quarks. 
We examined the simplest class of models where VLQs arise from a single representation 
under the SM gauge group. There are 7 possible representations where the VLQs have 
Yukawa interactions with the SM third generation quarks and Higgs doublets. The new 
complex Yukawa couplings could accommodate new sources of CP violation. Among them, 
we find that an irreducible CP phase shows up only for one representation of VLQ, which is 
a doublet under SU(2)l and carries hypercharge 1/3. For the other representations all the 
phases can be rotated away and are unphysical. 

We study the CP violating phenomenology in the doublet VLQ model. Since the VLQs 
are already constrained to be heavier than 800 GeV by the LHC, we integrate them out and 
study the effective theory, where CP violation manifests itself through a new right-handed 
charged current mediated by the SM IF-boson. We have calculated the CP violating Higgs 
interactions with SM gauge bosons, generated at one loop level involving both top and 
bottom quarks. This corresponds to a dimension 8 operator in the heavy top/bottom quark 
limit. As a consequence, only the hWW coupling is CP violating, while the hZZ, / 177 , 
hZ'y couplings are essentially CP conserving at this order. The strength of the CP violating 
hWW coupling is proportional to the quantity Iui^Or), where cll,r are the coefficients 
of left- and right-handed current Wtb interactions, respectively. At low energy, we fold 
the most relevant constraints on Im(a£a^) come from the electric dipole moments and the 
b —> sy decay rate and CP asymmetry, which are in complimentary to each other. The 
current experimental constraints require lm(a* L aR) < 0.01. They in turn imply that the 
coefficient of the CP violating hWW interaction, , cannot be larger than of order 10~ 5 , 
and, as we stress again, only in the hWW channel. We expect exciting interplays of various 
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experimental searches in the future to probe and distinguish new sources of CP violation 
near the electroweak scale. 
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